Glass Transition in a Two-Dimensional Electron System in Silicon in a Parallel 

Magnetic Field 



J. Jaroszynski/'lj Dragana Popovic, 1 and T. M. Klapwijk 2 

1 National High Magnetic Field Laboratory, Florida State University, Tallahassee, Florida 32310 
Department of Applied Physics, Delft University of Technology, 2628 CJ Delft, The Netherlands 

(Dated: February 2, 2008) 

Studies of low-frequency resistance noise show that the glassy freezing of the two-dimensional 
electron system (2DES) in Si in the vicinity of the metal-insulator transition (MIT) persists in 
parallel magnetic fields B of up to 9 T. At low B, both the glass transition density n g and n c , 
the critical density for the MIT, increase with B such that the width of the metallic glass phase 
(n c < n 3 < n g ) increases with B. At higher B, where the 2DES is spin polarized, n c and n g no 
longer depend on B. Our results demonstrate that charge, as opposed to spin, degrees of freedom 
are responsible for glassy ordering of the 2DES near the MIT. 

PACS numbers: 71.30.+h, 71.27. +a, 73.40.Qv 



The fascinating strong correlation physics exhibited by 
low-density two-dimensional (2D) electron and hole sys- 
tems Q remains the subject of intensive research. In the 
vicinity of the apparent metal- insulator transition (MIT), 
in particular, both electron-electron interactions and dis- 
order appear to be equally important. Their competition 
may lead to the emergence of many metastable states 
and the resulting glassy dynamics of electrons. Recent 
experiments 0, Q on a 2D electron system in Si have 
demonstrated such glassy behavior, lending support to 
the theoretical proposals that attempt to describe the 
2D MIT as the melting of a Coulomb 0,0, Wigner 0, 
or spin 7J glass. Even though several features of the 
data 0, 3 are consistent with the model ofglassy be- 
havior that occurs in the charge sector d, H, 0, it is 
still an open question whether charge or spin degrees of 
freedom are responsible for the observed glass transition. 
Since a sufficiently strong magnetic field is expected to 
destroy the spin glass order 0, , experimental stud- 
ies of glassy dynamics in parallel magnetic fields B ^l| 
should be able to distinguish between the proposed mod- 
els. Here we present such a study, which shows that the 
glass transition persists even in B such that the 2D sys- 
tem is spin polarized. These results demonstrate that 
charge, as opposed to spin, degrees of freedom are re- 
sponsible for glassy ordering of the 2DES near the MIT. 

Previous studies carried out at B — employed a com- 
bination of transport and low-frequency resistance noise 
measurements @, li| to probe the glassy behavior. The 
glass transition was manifested by a sudden and dramatic 
slowing down of the electron dynamics and by an abrupt 
change to the sort of statistics characteristic of compli- 
cated multistate systems, consistent with the hierarchical 
picture ofglassy dynamics. It was also found [2j that the 
glass transition occurs in the metallic phase, i.e. at an 
electron density n g > n c , where n c is the critical density 
for the MIT determined from the vanishing of activa- 
tion energy in the insulating regime |l2l The in- 
termediate metallic glass (MG) phase is of considerable 



width in strongly disordered samples {{n g — n c )/n c k, 0.5 
0), whereas in devices with low disorder n g is only a 
few percent higher than n c Q, in agreement with the- 
ory [9J. In this work, we investigate the glass transition 
in these same high peak mobility (low disorder) devices 
using noise spectroscopy in a parallel B. We find that 
all the signatures of the glass transition in parallel B are 
qualitatively the same as in B = 0, even when the elec- 
trons are spin polarized at high B. We construct a phase 
diagram in the (n s , B) plane (n s is the electron density), 
and find that the MG phase is broadened by a parallel B. 
The temperature (T) dependence of the conductivity a 
in the MG at low B is similar to the one first observed in 
the MG phase of highly disordered devices at B = , 
and it is consistent with theory ||. We show that, at 
low B, it also provides evidence for a quantum phase 
transition (QPT) at n c (B). 

The experiment was performed on n-channel Si metal- 
oxide-semiconductor field-effect transistors (MOSFETs) 
with the peak mobility /i «s 2.5 m 2 /Vs at 4.2 K, fab- 
ricated in a Hall bar geometry with Al gates, and ox- 
ide thickness d ox = 147 nm 14] . The resistance R was 
measured using a standard four-probe ac technique (typ- 
ically 2.7 Hz) in the Ohmic regime. A precision DC volt- 
age standard (EDC MV116J) was used to apply the gate 
voltage, which controls n s . We note that n s was always 
varied at a relatively high T«2K. Contact resistances 
and their influence on noise measurements were mini- 
mized by using a split-gate geometry, which allows one 
to maintain high n s (w 10 12 cm~ 2 ) in the contact region 
while allowing an independent control of n s of the 2D sys- 
tem under investigation in the central part of the sample 
(120 x 50 /im 2 ). The samples and the measurement tech- 
nique have been described in more detail elsewhere 

For each n s and B, R was measured as a function of 
time t at T = 0.24 K, although measurements at higher 
T were also performed at several selected B. At B = 0, 
the temperature coefficient of the time-averaged resistiv- 
ity d(p)/dT = at n* ss 9.7 x 10 10 cm- 2 , similar to what 
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FIG. 1: (a)AR(t)/{R), and (b) the corresponding power spec- 
tra Sn(f), at several B for n s = 11. 9x 10 10 cm -2 . In (a) traces 
are shifted for clarity. In (b) Sr(/) are averaged over octaves 
and multiplied by /, so that 1/f spectrum is horizontal on 
this scale. Solid lines are linear least-squares fits with the 
slopes a = 1.58, 1.24, 0.96, 0.58 (from top). 



was obtained on the previous cooldown of the same sam- 
ple 3] . At a fixed T and in the range of n s under investi- 
gation, (p) exhibits a dramatic increase with £?, followed 
by a much weaker dependence ( "saturation" ) at higher 
fields (B > 2 — 4 T). This large positive magnetoresis- 
tance at low B has been observed and studied extensively 
in many 2D systems [l| , including other samples from the 
same source u^i as ours. In the saturation region, it has 
been shown [15j that the 2DES is spin polarized. 

Figures nja) and^b) show the time series of the rela- 
tive changes in resistance AR(t) / (R) and the correspond- 
ing power spectra Sr(J), respectively, for a fixed n s and 
several B. It is obvious that B has a strong effect on both 
the amplitude and the character of the noise, as discussed 
in detail below. In order to compare the noise magnitudes 
under different conditions, the power Spt,(f = 1 mHz) is 
taken as the measure of noise. It is determined from the 
fits of the octave-averaged spectra to the form 1/ f a for 
1CT 4 < / < 0.07 Hz [solid lines in Fig. []]»]. In addi- 
tion, we have also analyzed the so-called second spectrum 
SaC/z)/)) which is the power spectrum of the fluctua- 
tions of Sii(f) with time 0|. S , 2(/2, /) provides a direct 
probe of correlations between fluctuators: it is white (in- 
dependent of f2) for uncorrelated, and S2 <x I//2 for 
interacting fluctuators 16] . At B = 0, the glass tran- 
sition in Si MOSFETs was manifested by a sudden and 
dramatic increase of Sr, a rapid rise of a from w 1 to 
a 1.8 2,3], and a change of the exponent (1 — (5) from a 
white (zero) to a nonwhite (nonzero) value [3j . We adopt 
similar criteria for the glass transition in B. 

Figures E^a), 13b), and|^c) present the i3-field depen- 
dences of Sr, a, and (1 — (3), respectively, for several n s . 
For each n s and B, ^(j^, /) was calculated for three 
different octaves /: 2-4 mHz, 4-8 mHz, and 8-16 mHz. 
In order to reduce the uncertainty in (1 — /3), the expo- 
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FIG. 2: (a) Sn(f = 1 mHz), (b) a, and (c) (1 - /3) vs. B 
for n s (10 10 cm" 2 ) shown on the plots; T = 0.24 K. Other 
data have been omitted for clarity. £/?(/) has been corrected 
for the white background noise. The arrow in (a) shows B g 
for n s (10 10 cm -2 ) = 11.9. The error bars on the right show 
maximum standard deviations of the data, which, for clarity, 
were plotted after performing a simple 3-point average. The 
origin of the fluctuations of Sr(B) for some n s at B > B g (e. 
g. for n s (10 10 cm -2 ) = 11.2) is not understood at this time. 



nent shown in Fig. Efc) represents the average (1 — j3) 
obtained from S2 in those three octaves. 

At the highest n s , in the metallic phase, "pure" 1/f 
noise (a ~ 1) resulting from uncorrelated [(1 — 0) w 0] 
fluctuators is observed in B = 0, Q . A decreasing Sr 
with B in Fig.^a) shows that parallel B suppresses such 
noise, the behavior that is already apparent from the raw 
data [e.g. bottom two traces in Fig. [Ja)]. In addition, 
the exponent a is reduced to ~ 0.5 [Fig. Hfb)]. The 
dependence of both Sr and a on B becomes weaker for 
B > 3—4 T. We note that noise in the metallicphase does 
not depend on n s , similar to the B = case 0, ■ Even 
though the origin of this noise is not known at the mo- 
ment, these data should provide a valuable contribution 
towards understanding the nature of the metallic phase. 
Pure 1/f noise in B = has been observed recently also 
in a 2D hole system in GaAs near the apparent MIT, on 
the metallic side 
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FIG. 3: (a) S R (f = 1 mHz), (b) a, and (c) (1 - 0) vs. n s for 
several B shown on the plot. The same symbols are used in 
all three panels. Dotted lines guide the eye. The maximum 
error bars are shown in the upper right. 

At lower n s , the noise at low B behaves as described 
above. However, for each n s , there is now a well-defined 
field B g [Fig. Eta)] where, after the initial decrease with 
B, Sr begins to increase dramatically, accompanied by a 
rapid rise of a and an increase of (1 — p) to nonwhite val- 
ues, indicative of the onset of strong correlations. This 
striking change of both the magnitude and character of 
noise within a narrow range of fields is obvious even from 
the raw data (e.g. middle two traces in Fig. In anal- 
ogy with the B — case, we identify B g as the field 
where glass transition occurs at a given n s . At even 
higher fields (B > 4 T), S R , a, and (1 - 0) no longer 
depend on B, consistent with the fact that the electrons 
are spin polarized. Nevertheless, Fig. El shows that the 
strong dependence of noise on n s is still present for a 
given B in this regime. The dependence of Sr, a, and 
(1 — pi) on n s is plotted explicitly in Fig. [3|for several B. 
The density where Sr begins to increase dramatically, 
accompanied by changes in a and (1 — /?), is identified as 
the glass transition density n g for a given B. 

The values of B g (n s ) and n g (B) determined in this 
way have been used to construct a phase diagram in the 
(n s ,B) plane, as shown in Fig. 0] The square symbols 
designate the boundary of the glassy phase, and clearly 
show an increase of n g with B, followed by a saturation 
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FIG. 4: T — phase diagram. The dashed lines guide the 
eye. Squares: boundary between the metallic and MG phases; 
other symbols: boundary between the MG and insulating 
glass phases. The data from Ref. [l3f have been shifted up 
by 0.85 x 10 10 cm -2 to make the n c (B = 0) values coincide. 



at higher fields. It is interesting to compare the behavior 
of n g (B) with n c (B). In Ref. [Rj, where samples almost 
identical to ours were used, n c (B) was determined based 
on both a vanishing activation energy and a vanishing 
nonlinearity of current-voltage characteristics when ex- 
trapolated from the insulating phase. We have used the 
activation energy method but the range of accessible T 
and n s in our ac measurements was smaller compared 
to the dc technique and dilution refrigerator T (down to 
30 mK) used in Ref. [13]. We find that the data that 
we have available in the insulating regime are best de- 
scribed by (a) oc exp[— (To/T)"] with n = 1/2, which 
corresponds to variable-range hopping with a Coulomb 
gap |l9| . The extrapolation of To(n s ) to zero, where 
only n s with To > 0.5 K were used, was employed to 
determine n c shown in Fig. ^ Strictly speaking, the 
rather limited range of data does not allow one to make 
an accurate distinction between different forms of acti- 
vated (a(T)). This experimental uncertainty is reflected 
in rather large error bars for n c shown in Fig. 0] Nev- 
ertheless, the agreement between our results and those 
obtained in Ref. 0] is remarkably good (Fig. 0}. 

Furthermore, it was possible for the first time to de- 
termine n c in a magnetic field by studying (<j(T)) on 
the metallic side of the MIT. In particular, in a rela- 
tively narrow range of n s for B = 1,2,3 T, the data 
are best described by the metallic power-law behavior 
(a(n„B,T)) = (a(n s ,B,T = 0)} + b(n s , B)T 15 [Fig. 
Eta)], similar to what was observed in the MG phase of 
highly disordered samples at B = |2| ■ The extrapolated 
T = conductivities go to zero at n c (B) in a power-law 
fashion (a(n s ,B,T = 0)) cx 5% with fi ~ 1.5 [see Fig. 
Etb); S n — n s /n c (B) — 1], in agreement with theoretical 
expectations near a QPT |2fJ. At B = 0, there is some 
evidence [2l| of similar behavior with fj, ~ 1 — 1.5, ob- 
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FIG. 5: (a) (a{T)) in the MG phase for n s (10 10 cm- 2 ) = 
11.9,11.6,11.3,11.2,11.0,10.9,10.7 from top; B = 2 T; 
n c (B = 2 T) = 10.67 x 10 10 cm- 2 . (b) (<x(T = 0)) oc $£. 



taincd from <j{T = 0) in the metallic phase of different 
Si MOSFETs. Considering that n c {B) for B = 1,2,3 T 
shown in Fig. g]have been obtained by approaching the 
MIT from the metallic side, the agreement between our 
results and those of Ref. [l3( is even more remarkable. 

The phase diagram in Fig. 0] makes it clear why the 
metallic T 3 / 2 temperature dependence of conductivity 
[Fig. Eta)] is observed in such a narrow range of n s : it 
is characteristic of the metallic glass phase, even in B of 
up to 3 T. At B — 0, the MG phase is very narrow and 
the T 3 / 2 correction is, therefore, difficult to observe in 
these samples, in contrast to highly disordered ones Q. 
A parallel B increases the width of the MG phase (Fig. 
gj. Inaddit ion, the range of T where T 3 / 2 holds increases 
with parallel B (not shown). The increase of n g and n c , 
and the broadening of the MG phase with B can be un- 
derstood to result from the suppression of screening by 
a parallel B which increases the effective disorder. 
This, in turn, favors glassiness, consistent with theoret- 
ical expectations 0- It is also interesting to note that 
the prefactor b{n s , B) of the T 3 / 2 correction [slopes in 
Fig. Efa)] does not seem to depend on n s , in contrast 
to the strong n s -dependence seen in highly disordered 
(low-mobility) samples at B = @. Further careful in- 



vestigation is required in order to determine whether this 
difference can be attributed to the effects of B or to the 
effects of disorder. Such a study, however, along with 
a detailed analysis of (a(T)) in the MG phase at high 
fields, is beyond the scope of this paper. 

In summary, the glass transition in high-mobility Si 
MOSFETs has been found to persist in high parallel B, 
where the 2DES is spin polarized. These results demon- 
strate that charge, as opposed to spin, degrees of freedom 
are responsible for glassy ordering of the 2DES near the 
MIT, consistent with recent theory 0,0,0. 
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